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INTRODUCTION 

In  certain  applications  It  Is  desirable  to  perform 
a  continuous  measurement  of  the  frequency  of  a  sine  wave 
when  the  sine  wave  Is  Immersed  In  a  random  noise  background. 

This  Is  a  common  application  of  systems  such  as  phase  lock 
loops,  as  well  as  other  forms  of  frequency  modulation 
detectors.  If  the  frequency  of  the  sine  wave  Is  very  stable 
over  a  long  period  of  time,  then  Its  frequency  can  be 
measured  with  arbitrary  accuracy.  However,  If  the 
frequency  fluctuates  In  some  random  manner  with  time 
and  the  goal  Is  to  track  these  frequency  fluctuations  In  the 
presence  of  noise,  then  an  optimum  system  can  be  derived  for 
estimating  the  frequency  versus  time  and  a  lower  bound  can 
be  placed  on  the  mean  squared  tracking  error.  This  Is  the 
topic  to  be  explored  In  this  report,  with  the  goal  of  providing 
insight  Into  the  design  and  performance  of  frequency  trackers. 
While  the  analytical  approach  used  Is  based  on  continuous 
Kalman  filtering  theory?  the  resulting  system  closely  resembles 
a  conventional  phase  lock  loop.  Thus  the  results  can  be  applied 
to  the  design  of  an  adaptive  phase  lock  loop  tracking  system. 

•The  notation  used  In  this  report  follows  that  in  Reference  (1) 
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ANALYTICAL  MODEL 

If  the  signal  to  be  tracked  is  a  sine  wave  with  instantaneous 
amplitude  A  and  radian  frequency  w,  it  may  be  represented  at 
any  instant  as  A  cos  8  where  0  ■  /in  dt,  The  instantaneous 
state  of  the  input  signal  process  may  be  represented  by  a  three- 
dimensional  vector 

0  phase  in  radians 

X]  ■  to  *  frequency  in  rad/sec  (1) 

A  amplitude  in  volts 

__  m 

The  frequency  parameter  u  is  the  state  variable  which  we 

wish  to  measure,  while  8  and  A  are  required  to  define  the  measure 

h(X)  -  h(  0  ,<*> ,  A)  =  A  cos0  (2) 

from  which  we  must  infer  the  state.  If  a  measurement  z  is  taken 

of  the  input  signal  plus  noise  process,  then  z  *  h(X)  +  v  where 

v  is  a  zero-mean  random  process  representing  the  contamination 

of  the  measurement  due  to  noise.  If  the  noise  is  a  white  noise 

process  with  a  (double  sided  frequency)  power  density  of  N 
2 

volts  /Hz,  a  noise  autocorrelation  matrix  R  can  be  defined  as 

v(t)  v(t+x)T  ■  [R]6 (t )  ■  N  6 ( t )  (3) 

where  R  becomes  a  lxl  matrix  equal  to  N. 


5 


NSWC/WOL/TR  75— 71* 


The  evolution  of  the  state  of  the  input  process  with  time 
may  be  described  by  the  set  of  differential  equations 

X]  -  [F]*X]  +  [Q]*w]  (4) 


where  X]  Is  the  first  derivative  of  the  state  X].,  the  matrix 


defines  the  dynamic  behavior  of  the  input  process,  w]  is  a 
random  white  noise  vector  which  produces  the  random 
fluctuations  in  the  input  signal,  and  [0]  is  a  matrix  which 
couples  this  random  process  to  the  system.  The  random  vector 
w]  is  described  by  its  mean  wQ]  and  a  covariance  matrix  [Q]. 
In  order  to  provide  random  fluctuations  in  both  frequency 
and  amplitude,  w]  is  defined  as  a  two-dimensional  vector  such 
that 


[w( t )  -  wQ3[w(t+T)  -  w0]T  -  CQ]6(t) 
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The  matrix  G  is  defined  as 


G  - 


'0 

a 

0 


0 

0 

6 


(8) 


With  these  definitions,  and  become  the  mean 
frequency  and  the  mean  amplitude  of  the  sinusoidal  process 
respectively,  while  Q  and  Q.  respectively  determine  the  sizes 

UJ  A 

of  independent  frequency  and  amplitude  perturbations  around 
these  means.  As  a  result  of  the  form  of  [P],  these  perturbations 
both  resemble  simple  low-pass  random  processes  which  may  be 
described  by  their  autocorrelation  functions 


♦  uw(t)  “  (aQw/2)  exp  (-a  |  t  | ) 


(9) 


and 


4>aa(t)  "  (BQfl/2)  exp (-8 1 t  | ) 


(10) 


The  rms  frequency  deviation  from  the  mean  u> o  is' thus 
AlQ^/2  radians/sec,  and  the  correlation  time  constant  (or 
typical  duration)  of  these  fluctuations  is  1/a  seconds? 
Similarly  the  rms  deviation  of  the  amplitude  from  the  mean 


Aq  is  volts  and  the  typical  time  associated  with  these 


*a  may  also  be  called  the  bandwidth  of  the  random  frequency 
fluctuations,  but  one  must  be  careful  to  distinguish  this 
modulation  "bandwidth"  from  the  rms  frequency  deviation  of 
the  signal. 
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fluctuations  Is  1/6  seconds.  While  this  model  Is  not 
Intended  to  represent  any  particular  system,  it  is 
sufficiently  general  to  represent  a  wide  range  of  sinusoidal 
processes  having  some  degree  of  both  amplitude  and  frequency 
random  modulation.  It  contains  sufficient  parameters  to 
characterize  both  a  magnitude  and  a  bandwidth  (or  typical 
period)  associated  with  each  type  of  fluctuation. 
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DERIVATION  OF  KALMAN  FILTER 

The  model  presented  in  the  previous  section  is  that  of 
a  continuous  linear  system  described  by  a  set  of  first 
order  differential  equations,  on  which  a  single  continuous 
measurement  is  made  where  the  measurement  is  a  nonlinear  function 
of  the  system  state.  The  theory  of  Kalman  filtering 
provides  a  means  of  estimating  the  state  of  this  system 
as  a  function  of  time,  using  the  noisy  input  measurements, 
such  that  the  mean  squared  error  in  the  estimate  is  minimized. 

A 

The  Kalman  filter  provides  a  continuous  estimate  X]  of  the  current 
state  X]  of  the  input  process  plus  a  covariance  matrix  [P] 
representing  the  estimated  uncertainty  in  the  estimate  of  X]. 

In  the  present  problem 


*  A  rn 

[x  -  x][x  -  xr 


ee 

P0o) 

TJ 

CD 

> 

00) 

P 

0)0) 

Po)A 

6A 

Po)A 

PAA 

(id 


contains  the  estimated  variance  of  the  estimate  of  each  state 
variable  plus  the  covariance  between  them.  Reference  (1)  gives 
the  form  of  the  continuous  Kalman  filter  as  a  system  obeying  the 
differential  equations 


X  «  F  X  +  G  WQ  +  P  HT  R"1  (z  -  h(X) ) 


(12) 


*  T  T  T  -1 

P  «  F  P  +  P  P1  +  G  Q  G1  -  P  H1  R  1  H  P 


(13) 
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The  H  matrix  in  this  set  of  equatiov  3  is  a  linearized  form 
of  the  measurement  function  and  consists  of  the  Jacobian  matrix 

[H]  -  dh/dX  (14) 

In  our  problem  this  becomes  a  1x3  matrix 


H  *  [3h/30,  3h/3w,  3h/3A] 


^  A  A 

*  [-A  sin0,  0,  cosG] 


(15) 


Thus,  given  an  initial  estimate  X(0)  of  the  state  and 
an  estimate  P(O)  of  the  initial  uncertainty,  these  matrix 

A 

equations  provide  a  way  of  computing  the  best  estimate  X 
as  a  function  of  time  as  well  as  showing  the  way  the 
uncertainty  matrix  evolves  with  time.  Note  that  the 
equation  for  the  behavior  of  P  depends  only  on  the  several 
matrices  defining  the  problem  and  not  at  all  on  the  measurements 
z  of  the  input.  The  differential  equations  for  the  state 

-■  A 

estimate  X  do  use  the  measurement  inputs  and  require  P  in 
their  solution.  Thus  the  general  form  for  the  Kalman  filter 

A 

consists  of  a  time-varying  filter  which  determines  the  X 
estimates  from  the  input  data,  where  the  coefficients  of  this 
filter  vary  with  the  estimated  error  matrix  P.  However 
in  problems  such  as  the  one  studied  here  the  system  model  and 
all  noise  matrices  are  stationary  in  time  so  that  a  steady-state 
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value  will  eventually  be  reached  for  the  P  matrix.  A  steady- 
state  form  of  the  filter  can  thus  be  found  by  setting  P 
equal  to  zero  In  equation  (13) ,  solving  for  Pt  and  then 
using  this  steady-state  P  matrix  In  equation  (12)  for  the 
state  estimation  filter. 


11 
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EXPANSION  OF  KALMAN  FILTER  EQUATIONS 

In  principle  the  freqeuncy  tracking  problem  is  solved 
at  this  point.  A  time-varying  Kalman  filter  may  be  implemented 
directly  by  repetitive  matrix  solution  of  equations  (12)  and 
(13)  using  the  matrices  as  defined  in  the  model  and  using  the 

A 

residuals  (z-h(X))  of  the  input  measurements  as  inputs  to 
the  filter.  In  large  Kalman  filters  (having  a  large  number 
of  state  and/or  measurement  variables)  this  direct  matrix 
approach  is  about  the  only  method  available  because  analytical 
expansion  of  the  deceptively  simple  looking  matrix  equations 
is  impractical.  In  this  rather  small  problem,  however,  it 
is  possible  to  expand  the  indicated  matrix  operations 
analytically  and  gain  much  insight  into  operations 
actually  being  performed  by  the  filter. 

The  first  necessary  step  is  to  expand  the  matrix  equation  (13) 

t 

for  P  to  form 


P  P  +  P  FT  +  G  Q  0T  -  P  HT  R"1  H  P 


P6(i) 

P 

(1)0) 

PtoA 

- 1 

*TJ 

<D 

e 

”ap0(l) 

"6PeA 

- 

“aP(jJU 

”aPwA 

+ 

Pww 

”aP(iKU 

"BPuA 

+ 

.  -eP6A 

-BPUA 

-6Paa  _ 

_P“A 

-°PUA 

-bpaa 

% 

0 

0 


-APee  sin6  +  p0A  cose 

-AP0q  sine  +  PQA  cos6 

(1/N) 

•APeu  8lne  +  pwA  0080 

• 

-AP0U  sln0  +  PuA  oosO 

! -AP0a  sin9  +  P^A  cos0 

-AP0A  8in0  *  PAA  0080 

m 

12 


NSWC/WOL/TR  75-74 


From  this  matrix  equation  we  can  extract  the  six  scalar 
differential  equations  for  the  components  of  P 

P00  -  2PQu  -  (l/N)(-APe0sine+P0Acose)2  (17a) 

+  -  (l/NH-AP^slne+P^oose)2  (17b) 

PAA  -  -26?  ^  +  62Qa  -  (l/N)(-APeAsine+PAAcose)2  (17c) 

P9u  '  ?m  -  “P9U  -  a/NH-APeeSinew^coseH-AP^slne+P^cosO)  (17a) 

P8A  ■  PuA  -  6P0A  -  (l/N)(-AP69sln6+PeAco8e)(-APeA8lne+PAAcose)  (17e) 

PwA  *  “(«+B)pwA  -  (l/N)(-AP0ajsine+PajAcos8)(-AP6Asin0+PAAcos0)  (17f) 

Solving  these  six  differential  equations  continuously 

and  exactly  is  a  rather  awesome  task;  fortunately  some 

simplifications  are  possible.  First,  tnere  is  an  implicit 

time  dependence  in  the  coefficients  of  these  equations  since 

sin0  and  cos0  both  vary  at  the  input  signal  frequency. 

This  time  dependence  is  induced  by  the  dependence  of  the 

measurement  matrix  H  on  the  signal  phase  0.  If  we  assume  that 

fluctuations  in  the  P  matrix  (that  is,  variations  in 

uncertainty  of  the  state  estimate)  occurring  within  a  single 

period  of  the  input  sinusoid  are  unimportant,  then  time-average 

values  of  P  may  be  used  in  determining  the  dynamics  of  P, 

where  the  time  average  is  taken  over  a  period  of  the  input 

2  2 

sinusoid.  Using  <cos  0>  *  <sin  0>  ■  1/2  and  <sin0  cos0>  *  0, 
the  time-average  form  of  equations  (17)  becomes 
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<%e>  -  2Peu  -  a/2NHA2P^6  ♦  p^a>  (18a) 
<*««>  *  -2aPuu  +  a\  '  C1/2N) (A2PgM  +  P2a)  (18b) 
<PAA>  =  -2flPAA  +  ®\  -  <1/2N)(A2P2a  ♦  P2a>  (18c) 

<N«>  -  p»  -  °Pe„  -  (l/2H)(A2PeePeu  ♦  P6A  PuA)  (188) 
‘'•A"  '  P»A  -  ®P6A  *  (l/2NXA2pe9p6A  ♦  P6A  PAA)  (180 
<PuA>  .  -(a+S)PuA  -  (1/2N)(A2P6uPm  ♦  PuA  PAA)  (I8f) 


A  second  simplification  is  now  possible  by  recognizing 
that  equations  (l8e)  and  (l8f)  are  linear  in  P0A  and  P^A  and 
that  there  are  no  other  terms  driving  them  as  inputs.  Their 
solutions  are  decaying  exponentials  so  that  both  of  these 
covariances  decay  to  zero  if  for  any  reason  they  are  initially 
given  non-zero  values.  This  is  physically  reasonable  since 
there  is  no  coupling  mechanism  in  the  model  between  amplitude 
and  phase  or  frequency.  Thus  one  would  expect  their  estimation 
errors  to  be  uncorrelated.  Note  that  this  reasoning  does 
not  apply  to  PQw  since  errors  in  frequency  naturally  evolve 
into  errors  in  phase.  However  the  fact  that  P0A  and  P^A  can 
be  assumed  equal  to  zero  reduces  the  number  of  equations  to 
four  and  simplifies  those  remaining  to 


• 

<pee>  * 

2P0 

00) 

-  (A2/2N)P2Q 

(19a) 

<p  >  - 
0)0) 

a\ 

-  2oP« 

-  (A2/2H)p|B 

(19b) 

• 

<P  >  = 
rAA 

s2«a 

-  26Paa 

-  <i/2N>paa 

(19c) 

• 

<PQ  >  » 
00) 

P 

0)0) 

'  aPeM  - 

(A2/2N)pe0Peu 

(19d) 
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Having  found  the  differential  equations  representing 
the  behavior  of  the  uncertainty  matrix  P  with  time,  we  can 
proceed  in  a  similar  fashion  with  the  state  estimation 
equation,  equation  (12).  Substituting  and  multiplying  the 
indicated  matrixes  in  equation  (12)  gives 


X«FX  +  Gw0  +  PHT  R_1(z  -  h(X) ) 


• 

0 

l 

0 

A  m 

e 

’o 

0  * 

0J 

»  « 
-APggSine 

s 

0 

-a 

0 

A 

•  0) 

A 

+ 

a 

0 

0 

• 

A„ 

+  (1/N) 

-APeuslne 

0 

- 

0 

-e 

- 

A 

m 

0 

e 

0 

.  PAAoose  , 

A  A 

•(z-A  cose) 


(20) 

where  we  have  taken  advantage  of  the  fact  that  PgA  and  PyA 
are  equal  to  zero.  Writing  tliese  out  as  individual  differential 
equations  gives 


•  A  A  A 

6  *  u  -  (A/N)P0gSine  (z-A  cos0) 

A 

•  A  A  A 


(i)  *  o(u0-w)  -  (A/N)Pj  ^sinO 

A 

A  =  B(A0-A)  +  (1/N)Paacos0 


(z-A  cos6) 

A  A 

(z-A  cos0) 


(21a) 

(21b) 

(21c) 


Just  as  in  the  equations  for  P,  these  equations  contain  some 
terms  which  vary  sinusoidally  at  the  signal  frequency  or 
higher.  Assuming  that  these  short  term  variations  are  meaningless 
as  far  as  estimation  of  the  state  is  concerned,  we  can  again  take 
time  averages  of  these  equations  over  the  sinusoidal  signal 
period  to  obtain 
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a 

•  A  A 

<e>  »  U)  -  (A/N)Pgg  <z  sin0> 

A 

•  A  A 

«d>  »  o(w0-w)  -  (A/N)  Pqw<z  sin0> 

A 

<A>  -  0(Ao-A)  +  (1/N)  PAAC<z  cos0>  -  (A/2)] 


Note  that  once  the  equations  are  written  in  this  form  the 

A 

concept  of  forming  the  residual  (z-h(X))  and  using  it  to 
update  the  estimate  vanishes  and  is  replaced  by  the  idea  of 

A  A 

using  the  mean  products  of  z  with  sin0  and  cosG  as  the  basic 
information  inputs  to  the  filter. 

Equations  (19)  and  (22) thus  become  the  equations  which  must 
be  solved  to  implement  the  time-variable  Kalman  filter,  where 
advantage  has  been  taken  of  the  lack  of  correlation 
between  amplitude  and  frequency  estimation  errors  and  of 
smoothing  of  the  equations  over  the  period  of  the  input 
signal.  Assuming  some  method  is  available  for  solving  the 
nonlinear  equations  (19)  to  obtain  P,  the  state  estimation 
equations,  (22),  can  be  implemented  by  using  a  system 
such  as  that  shown  in  Figure  1.  Each  integrator  output 
represents  one  of  the  state  variables,  while  its  input  is  the 
derivative  of  that  state  variable  formed  according  to  equations  (22) . 

Several  interesting  things  can  be  learned  from  the  form  of 
the  system  shown  in  Figure  1.  First  is  that,  while  the 
subsystem  for  estimating  amplitude  is  nearly  independent  of  that 
for  estimating  frequency  and  phase,  each  depends  on  the  other 
for  some  piece  of  information.  An  estimate  of  A  is  used  to 
establish  gain  coefficients  in  the  frequency  subsystem,  while  the 

A 

estimate  of  0  is  required  to  generate  the  product  <z  cos0> 
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(22a) 

(22b) 

(22c) 


FIG.  1  STATE  ESTIMATION  PORTION  OF  KALMAN  FILTER 
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used  as  the  Input  to  the  amplitude  subsystem.  Values  of 
Pe0,  P0w,  and  PAA  are  also  used  by  the  state  estimation  system 
to  establish  gain  coefficients.  Both  optimal  subsystems  have 
forms  which  are  commonly  used  in  systems  of  this  type.  The 
subsystem  for  estimating  A  is  really  a  simple  low-pass  filter 
(or  exponential  averager)  with  a  3  db  bandwidth  of  0  +  PAA/2N 
radians/second  (or  a  time  constant  equal  to  the  Inverse  of  this 
bandwidth).  It  is  driven  by  two  inputs,  one  of  which  is  the 
presumably  known  average  signal  level  AQ  and  the  other  is  the 

A 

measured  "in  phase"  component  of  the  input  process  <z  cos0>. 

In  most  systems  encountered  in  practice,  AQ  is  not  known  and 

A 

the  entire  estimate  of  A  must  be  derived  from  the  <z  cos0> 
input . 

The  subsystem  for  estimating  frequency  and  phase  is  also 
easily  recognized  as  a  second-order  phase  lock  loop.  The 
second  Integrator,  whose  input  is  0,  usually  consists  of  a 
voltage  controlled  oscillator  (VCO)  which  directly 

A  A 

produces  the  sin0  and  cos0  outputs.  Except  for  the  effect 
of  the  small  feedback  loop  around  the  first  Integrator,  the 
control  function  for  the  VCO  is  of  the  "proportional  plus 

A 

integral"  type  where  the  output  <z  sin0>  of  the  phase  detector 
is  applied  to  both  the  first  and  second  integrator  inputs. 

The  primary  effect  of  the  feedback  around  the  first  integrator  is  to 
limit  its  low-frequency  gain  to  remain  constant  below  a 
radians/second  rather  than  continuing  on  to  infinity  at  zero 
frequency.  The  mean  frequency  input  w0  also  serves  to  return  the 
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frequency  estimate  to  this  value  in  the  absence  of  other 
information.  Again  this  use  of  a  known  mean  frequency  is  seldom 
encountered  in  phase  lock  loop  systems.  However,  it  is  a  useful 
concept  whenever  any  a  priori  information  about  the  center 
frequency  of  the  sinusoidal  signal  is  available  and  provides 
a  means  of  "tuning"  the  detector  prior  to  lock  on,  or  of 
prompting  it  to  expected  changes  in  average  frequency  such 
as  from  doppler  effects  between  transmitter  and  receiver. 
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STEADY  STATE  FORM  OF  FILTER 

When  a  Kalman  filter  Is  used  In  problems  where  all  defining 
matrices  are  stationary  In  time  (or  where  short-term  averages 
are  stationary  as  In  this  problem),  a  steady  state  Is  eventually 
reached  where  P  becomes  constant  and  the  gain  coefficients  In  the 
state  estimation  part  of  the  filter  also  become  constant. 

Unless  the  behavior  of  the  system  before  the  steady  state  Is 
reached  Is  of  particular  Importance,  It  is  usually  satisfactory 
to  build  a  fixed  parameter  version  of  the  state  estimation 
filter  using  the  steady  state  values  of  P  and  to  avoid  completely 
the  necessity  of  Implementing  the  differential  equations  for  P. 
The  steady  state  values  of  P  also  indicate  the  mean  square 
errors  associated  with  the  fixed  parameter  estimation  filter 
in  the  steady  state  and  thus  serve  as  a  performance  measure. 

Finding  the  steady  state  values  of  P  requires  setting  the 
equations  (19)  for  P  equal  to  zero  and  solving  the  resulting 
simultaneous  nonlinear  equations.  Setting  each  of  equations  (19) 
to  zero  and  rearranging  slightly  yields 


2  P8ti> 

-  (A2/2N)  P2e 

(23a) 

'  2«pu0)  +  PL 

(23b) 

e\ 

*  28PAA  +  (1/2N>  PAA 

(23c) 

P 

Ufa) 

■  «peu  +  (A2/2N)peep9ll) 

(23d) 

Equation  (23c)  contains  only  PAA  and  can  be  solved  by 
itself  as  a  quadratic  in  PAA  to  yield 
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PAA  *  -2B  N  +  A  62N2  +  262NQa 


-  1] 


(24) 


where  the  second  root  has  been  dropped  as  clearly  extraneous. 
The  remaining  three  equations  can  be  solved  simultaneously  by 
using  (23a)  to  substitute  for  Pfi  in  (23b)  and  (23d)  to  give 


a\  ■  20  P«o  +  (A2/2N)3p20/4 


(25a) 


-  a(A2/2H)P2e/2  +  (A2/2N)2P39/2 


(25b) 


Then  substituting  (25b)  into  (25a)  gives 


*\  -  “2p  pe0  +  ap2pee  +  p3pee/J* 


ppee(a  +  pPee//2^ 


(26) 


where 


p  =  AV2N 


(27) 


is  a  signal  to  noise  ratio  (strictly  signal  power  to  noise 
power  density,  with  the  dimensions  of  Hertz)  parameter. 
Multiplying  by  p  and  taking  the  square  root  on  both  sides 
converts  equation  (26)  into  a  quadratic  of  the  form 


21 


NSWC/WOL/TR  75-74 


(pPee)2/2  +  a(pPee)  +  a/5^  -  0 

for  which  the  solutions  are 

PQe  ■  [-<*  +  «/a2  +  2a/pS^ ]/p 

-  (a/p)[/{  +  2  /pQ^/c^  -  1] 

where  only  the  single  real  positive  solution  has  been 
Equation  (25b)  may  now  be  used  to  find  Pwu  as 

Pu>ta>  -  op(o/p)2[  XT  2  /p^/o'  -  l]2/2 
+  p2(o/p)3[ /TT/ZfyT  -  l]3/2 

■  (a3/2p)/^+^~ /FSjJ'/a'  [/^~+_2_" -  l]2 

and  equation  (23a)  gives  PQw  as 

P0a)  -  (p/2)(a/p)2[/l  +  2  /pQ^/a  "  -  l]2 

■  (a2/2p)[*4  +  2  /piJ^Va"'-  l]2 


(28) 

(29) 

retained. 


(30) 

(3D 


Equations  (24),  (29),  (30),  and  (31)  thus  represent  the  steady 
state  solution  for  the  uncertainty  matrix. 
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These  equations  for  the  components  of  P  do  not  provide 
much  physical  insight  when  written  in  terms  of  Q.  and  Q  . 

tx  U) 

They  reduce  to  their  simplest  forms  when  some  new  variables 
are  introduced,  representing  a  signal  to  noise  ratio  and 
a  modulation  index  for  each  of  the  modulation  components. 

The  rms  frequency  deviation  has  been  shown  to  be  D  *  /aQ^/2  ,v 
and  a  frequency  modulation  index  can  be  defined  of  the  form 


Iw  -  D/a  -  (32) 

This  definition  of  the  modulation  index  is  consistent  with  the 

usual  (frequency  deviation)/ (modulating  frequency)  used 

in  describing  FM  systems  since  a  represents  a 

bandwidth  of  the  modulating  process.  The  modulation  index 

can  be  loosely  identified  with  an  rms  phase  deviation,  although 

the  phase  is  strictly  a  random  walk  process  in  this  problem 

with  no  finite  variance. 

Similarly  an  rms  modulation  percentage  can  be  defined  for 
the  amplitude  modulation  component  by  dividing  the  rms  amplitude 
fluctuation  /{JQA/2'  by  the  mean  amplitude  of  the  signal.  This 
gives  a  parameter 

IA  -  /BQ^7T/A  (33) 
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The  signal  to  noise  density  ratio  p  can  be 
nondimensionalized  for  each  modulation  component  by  dividing 
by  the  appropriate  modulating  process  bandwidth.  This  gives 

pu  -  p/a  -  A2/2aN  (3*0 

and 

PA  -  p/6  -  A2/26N  (35) 

Physically  these  may  be  considered  to  be  the  non-dimensional 
signal  to  noise  ratios  measured  over  the  "natural"  bandwidths 
of  the  processes  being  estimated. 

In  terms  of  these  new  variables,  the  equations  for  the 
P  matrix  components  may  be  rewritten 

Paa/A2  -  [A  *  2  IAPA  '  -  1]/PA  (36) 

pee  ■  CA  ~  2  -  U/ea  (37) 

(o2/2)/l  +  2  +  2  V5S?  -i:2/Pu  (38) 

p9u  ■  <«/2>[/l  +  2Iw^Pj'  -  1]2/PU  (39) 

If  the  steady-state  values  of  the  P  matrix  as  derived  in 
equations  (36)  through  (39)  are  inserted  as  the  appropriate 
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«  f 


coefficients  in  the  system  shown  in  Figure  1,  a  steady-state 
form  of  the  Kalman  filter  can  be  derived  as  shown  in  Figure  2. 
Figure  2a  shows  the  amplitude  estimation  portion  of  the 
system,  and  Figure  2b  shows  the  frequency  tracking  portion. 

The  form  shown  for  the  amplitude  estimator  identifies  it  as  a 
simple  low-pass  exponential  averager  with  a  time  constant  ta 
given  by 


4  ^ +  »**"*] 


(HO) 


This  averager  has  two  inputs;  one  is  the  a  priori  average  signal 

A 

level  Aq  and  the  other  is  the  measured  quantity  <2z  cos0>. 

The  weight  attached  to  the  measured  data  input  is 


-  [^7 


2iA2pA 


-  1]/ 


/T 


2iA2pA 


(Hi) 


while  the  weight  attached  to  the  a  priori  information  is 
1  -  Ga«  ta  and  Ga  are  plotted  in  Figure  3  as  a  function  of. 

PA  for  several  values  of  the  modulation  index  IA>.  Notice 
that  if  either  the  modulation  index  or  the  signal  to  noise  ratio 
is  small,  the  time  constant  approaches  the  natural  time 
constant  1/0  of  the  input  fluctuations  and  most  of  the  weight 
is  attached  to  the  a  priori  estimate.  In  the  limiting  case  the 
output  simply  becomes  equal  to  Aq  and  the  measured  amplitudes 
are  ignored.  For  large  values  of  IA  and  pA  most  of  the  input 

A 

weight  goes  to  the  measured  amplitude  <2z  cos6>  and  the  time 

constant  of  the  averager  becomes  shorter.  In  the  limit  this  time 

constant  approaches  the  asymptotic  form 
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FIG.  2A  AMPLITUDE  ESTIMATION  CIRCUIT 


FIG.  2B  PHASE  AND  FREQUENCY  ESTIMATION  CIRCUIT 


FIG.  2  STEADY  STATE  FORM  OF  KALMAN  FILTER 


AVERAGE  TIME  CONSTANT  (r.) 
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FIG.  3A  VARIATION  OF  TIME  CONSTANT 


NORMALIZED  SIGNAL  TO  NOISE  RATIO  (pA>  IN  DECIBELS 


FIG.  3B  VARIATION  OF  INPUT  DATA  WEIGHT 
FIG.  3  VARIATION  OF  AVERAGER  PARAMETERS  WITH  S/N  AND  MODULATION  INDEX 
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TA  -  l/IAB/5p^  for  IA?pA  »  1 


(42) 


The  frequency  estimator  shown  in  Figure  2b  forms  a  phase 
lock  loop  circuit  with  a  local  feedback  loop  around  the  first 
integrator.  This  local  feedback  loop  causes  the  first  integrator 
to  behave  like  a  low-pass  filter  rather  than  a  full  integrator. 

A  transfer  function  can  be  derived  for  the  phase  lock  loop 

A 

system  by  linearizing  the  output  <2(z/A)  sin0>  of  the  phase 

A 

detector,  such  that  this  output  is  approximately  6-0 
for  small  phase  errors.  The  transfer  function  then  becomes 


T(S)  .  .  a(K-i)«s  +  (Ka-1)«2 

0 (s )  2s2  +  2Kos  +  (K  -l)cr 


(43) 


where 


K  - 


(44) 


The  poles  of  this  second  order  system  are  located  at 

spole  *  (a/2 ) (-K  +  j A2  -  2)  (45) 

which  form  a  complex  pair  whenever  K  >  /P  .  It  is  convenient 
to  discuss  the  "bandwidth”  of  the  phase  lock  loop  system  which 
may  be  defined  as  the  radius  to  these  poles.  This  bandwidth  is 
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BW  -  a  /{Y?-l)/2  -  aXjSg*  (46) 

and  is  plotted  in  Figure  4  for  the  region  where  the  complex  poles 
occur.  It  may  be  seen  from  equation  (45)  that  for  large  K  the 
real  and  imaginary  parts  of  the  pole  location  are  equal  (that  is, 
the  system  is  about  70£  critically  damped)  and  that  the  damping 
increases  for  smaller  K  until  critical  damping  is  reached  for 
K  «  /T. 

As  the  parameter  K  becomes  smaller  than  ST*  the  inputs 
from  the  phase  detector  become  progressively  weaker,  until  for 
K  very  near  unity  all  that  is  left  is  the  local  feedback  loop 
around  the  <*>  integrator  which  causes  the  frequency  estimate  to 
settle  at  u>0.  This  causes  the  system  to  degrade  gracefully 
to  the  a  priori  estimate  of  the  average  frequency  when  there  is 
insufficient  information  from  the  input  to  permit  a  better 
estimate. 
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STEADY  STATE  PERFORMANCE 

The  steady  state  values  of  the  P  matrix  components  given  in 
equations  (36),  (37),  and  (38)  also  indicate  the  mean  squared 
errors  in  the  amplitude,  phase,  and  frequency  estimates  in  the 
steady  state.  Figures  5,  6,  and  7  represent  the  rms  errors  expected 
in  each  of  the  estimated  quantities,  in  terms  of  the  appropriate 
signal  to  noise  ratio  and  modulation  index.  Figure  5  shows 
the  relative  uncertainty  oA/ A  of  the  amplitude  estimate  versus 
PA  and  IA,  where  the  signal  to  noise  ratio  pA  is  expressed  in 
decibel  form.  Note  that  in  the  limit  of  low  pA  the  relative 
RMS  error  is  just  equal  to  the  modulation  index  IA.  This  is 
consistent  with  the  previously  noted  fact  that  the  estimate 
in  this  case  is  just  equal  to  AQ  so  that  the  entire  modulation 
shows  up  as  an  error  output.  As  the  signal  to  noise  ratio 
increases  the  filter  begins  to  track  the  amplitude  fluctuations 
so  that  the  rms  error  begins  to  decrease.  For  a  sufficiently 
large  signal  to  noise  ratio,  the  amplitude  estimation  error 
asymptotically  approaches  the  form 

oA/A  -  ^aa/A2'*  (Ia)1/2(2/pa)1/1<  for  pA  »  1/IA2  (47) 

This  somewhat  strange  dependence  on  IA  and  pA  is  due  to  the  fact 
that  the  time  constant  of  the  estimation  filter  is  also  a  function 
of  these  two  parameters.  Since  the  ratio  aA/A  has  the  same 
significance  at  the  output  of  the  amplitude  estimator  as  the 
modulation  index  IA  at  the  input,  an  improvement  factor  equal  to 
the  reduction  in  uncertainty  of  the  signal  amplitude  may  be  defined 
and  written  in  the  form 
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NORMALIZED  RMS  FREQUENCY  ERROR  (aw  /a ) 
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NORMALIZED  SIGNAL  TO  NOISE  RATIO  (f  w)  IN  DECIBELS 


FIG.  7  FREQUENCY  TRACKING  ERROR  VS.  S/Kl  AND  MODULATION  INDEX 
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(oa/a)/ia  ~  (2/iA2pA)1/1‘  -  ( 2e/iA2p)1/J 4 

for  IA2pA  >  >  1  (4R) 

This  factor  represents  the  amount  by  which  the  rms  uncertainty 
in  the  amplitude  may  be  reduced  by  using  the  Kalman  filter 
estimator,  relative  to  just  using  the  a  priori  estimate  AQ. 

Figure  6  shows  the  rms  phase  error  of  the  phase  and 
frequency  estimation  part  of  the  system.  The  signal  to  noise 
ratio  p,,  is  expressed  in  decibels  and  the  modulation  index  I 

OJ  U) 

appears  as  a  parameter.  Note  that  the  rms  phase  error  is  not  bounded 
at  low  signal  to  noise  ratios  but  continues  to  rise  indefinitely. 

This  is  also  a  reasonable  result,  since  the  system  shown  in 
Figure  2b  really  gives  up  trying  to  track  phase  at  low  signal  to 
noise  ratios  and  the  phase  is  a  random  walk  process  with  infinite 
steady  state  variance.  Once  the  rms  phase  error  exceeds  some  upper 
limit  it  is  inappropriate  to  discuss  the  estimation  system  as  a 
phase  lock  loop,  since  the  estimator  is  not  truly  "locked"  to  the 
input  signal  and  the  linearity  assumption  used  in  deriving 
equation  (43)  ceases  to  hold.  If  we  arbitrarily  choose  an  rms 
phase  error  of  one  radian  as  the  upper  limit  for  which  we  will 
consider  the  phase  lock  loop  "locked",  then  a  relationship 
of  the  form 


<  (Pu/2)3/2  +  <P«/2>1/2 


(49) 
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can  be  derived  as  a  criterion  for  lockon.  For  combinations 
of  modulation  index  and  signal  to  noise  ratio  violating  this 
requirement  the  rms  phase  error  exceeds  one  radian  and  the 
main  function  of  the  phase  detector  output  is  to  "pull  the 
frequency  estimate  toward  the  true  signal  frequency  without 
actually  following  it.  For  a  sufficiently  large  signal 
to  noise  ratio  the  rms  phase  error  approaches  an  asymptotic  form 

°e  =  <V1/,,(2/pu)3/8  for  po>  ”  1/8I“2  (50) 

The  rms  error  in  the  frequency  estimate  is  shown  in 
Figure  7  as  a  function  of  the  signal  to  noise  ratio  pu 
(expressed  in  decibels)  and  the  modulation  index  1^.  The  rms 
frequency  error  plotted  along  the  vertical  axis  of  Figure  7 
is  normalized  by  dividing  by  the  "natural"  fluctuation 
frequency  oi  to  make  it  dimensionless.  This  normalized  rms  frequency 
estimation  error  can  be  seen  to  approach  Iu  as  a  limit  for  small 
signal  to  noise  ratios  p  .  Since  lw  is  defined  as  the  rms 
frequency  deviation  D  divided  by  a,  this  means  that  in  the  limit 
of  low  signal  to  noise  ratio  the  rms  frequency  error  is  Just 
equal  to  the  rms  deviation  of  the  input  frequency  from  the  mean 
a j  ,  Again  this  result  is  physically  reasonable*  since  we  have  seen 
that  the  phase  lock  loop  degenerates  into  a  simple  circuit  whose 
output  is  simply  wQ  for  this  case.  In  the  limit  of  high  signal  to 
noise  ratios,  the  rms  frequency  estimation  error  approaches  the 
asymptotic  form 
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%  *  ^  =  C2)1/2(o.)1(I03)3/4<2/pm)1/8  for  pu  »  1/8IU2  (51) 

It  is  instructive  to  re-introduce  the  rms  input  frequency 
deviation  D  in  equation  (51)  so  that  the  rms  error  can  be 
expressed  in  terms  of  an  improvement  factor 

ow/D  ~  (32/Iw2pw)1/8  ■  (32ot3/D2p)1/8  for  Sp^2  »  1  (52) 

where  the  improvement  factor  is  the  ratio  of  the  rms  estimation 
error  to  the  input  frequency  deviation.  The  improvement 
factor  thus  represents  the  reduction  in  signal  frequency 
uncertainty  possible  by  using  the  Kalman  filter  tracker 
to  estimate  the  signal  frequency,  as  compared  to  using  only  the 
a  priori  estimate  w0  for  the  signal  frequency  and  accepting 
the  rms  deviation  D  as  the  uncertainty. 

It  is  interesting  that  the  improvement  factor  only  goes  as  the 
eighth  root  of  the  signal  to  noise  ratio,  so  that  increased  signal 
to  noise  ratio  improves  the  tracking  accuracy  only  very  slowly. 

This  result  is  partly  due  to  the  model  used  here  for  frequency 
fluctuations,  where  the  random  modulation  process  has  significant 
power  at  high  frequencies.  Thus  even  with  the  large  tracking 
bandwidths  occurring  at  high  signal  to  noise  ratios, 
frequency  fluctuations  of  still  higher  frequency  occur  and  are 
not  tracked.  A  signal  model  assuming  smoother  behavior  of  the 
frequency  fluctuations  would  not  necessarily  give  the  same  result. 
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The  shaded  area  of  Figure  7  represents  the  region  where  the 
rms  phase  error  exceeds  one  radian  and  the  loop  is  not  locked 
in  phase.  Since  the  linearized  Kalman  filter  derivation 
assumed  small  phase  errors,  this  portion  of  the  result  violates 
these  assumptions  and  is  not  necessarily  correct.  Some  tendency 
to  track  the  input  frequency  fluctuations  is  to  be  expected 
due  to  the  "frequency  pulling"  effect  in  an  unlocked  phase  lock 
loop,  but  it  is  not  likely  to  be  as  good  as  that  predicted  by  a 
linearized  analysis.  Thus  the  actual  performance  can  be  expected 
to  follow  the  solid  curves  in  the  phase  locked  region  but  to  degrade 
rapidly  toward  the  no-signal  performance  of  o^/ a  «  I  in  going 
through  the  marginal  phase  lock  region  toward  the  unlocked  region. 
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RECEIVER  LOCKON  THRESHOLD 

As  discussed  in  the  previous  section,  for  sufficiently 
low  signal  to  noise  ratios  the  rms  phase  errors  are  so  large 
that  the  linearity  assumptions  used  in  deriving  the  receiver 
are  violated.  Since  the  receiver  operation  clearly  relies 
on  a  phase  lockon,  it  will  fail  to  operate  properly  once  the 
phase  errors  become  excessively  large.  This  is  evidenced  by 
failure  to  track  the  signal  frequency  and  occasional  "cycle 
slips"  between  the  input  signal  and  the  receiver  oscillator. 

The  choice  of  the  maximum  permissible  phase  error  is 
somewhat  arbitrary  since  tracking  failures  simply  occur  with 
increasing  frequency  as  the  rms  phase  error  increases.  However, 
a  good  case  can  be  made  for  choosing  oQ  *  one  radian  as  the 
lockon  criterion.  Recovery  from  phase  errors  as  great  as 
two  radians  is  reasonably  certain,  but  a  phase  error  of  tt 
radians  generally  causes  a  cycle  slip.  Since  errors  of  2d 
occur  fairly  frequently  in  a  gaussian  process  while  3 o  errors 
are  comparatively  rare,  the  performance  of  a  phase  lock  receiver 
can  be  expected  to  degrade  rapidly  once  Oq  exceeds  one  radian. 

Equation  (49)  gives  the  relationship  between  Iu  and  Py 

for  an  rms  phase  error  below  one  radian.  This  can  be  converted 
to  more  physically  meaningful  terms  by  using  equations  (32)  and 
(34)  to  give 

» •  •  [®!/!  •  (*) ,/!] 
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In  practical  applications  it  is  useful  to  express  the 

signal  to  noise  ratio  (SNR)  as  the  signal  power  divided  by 

the  noise  power  per  Hertz  of  a  one-sided  (positive  frequencies 

only)  spectrum.  Since  p  was  defined  in  terms  of  the  double- 

sided  noise  density  N,  the  one-sided  noise  power  per  Hertz 

2 

is  2N  while  the  signal  power  is  A  /2,  so  we  have 


SNR  « 


a2/2 

2N 


P 

2 


(5*0 


Consequently  for  lockon  with  crQ  <  one  radian  we  have 


D  <  a  [(SNR/a)3/2  +  (SNR/a)1/2 ]  -  (a.SNR)1/2(l  +  SNR/a)  (55) 

This  result  is  plotted  in  Figure  8,  showing  the  minimum 
required  SNR  for  tracking  as  a  function  of  the  rms  frequency 
deviation  D  (expressed  in  Hertz  rather  than  radians/second)  for 
five  values  of  modulation  bandwidth  a,  For  a  given  value  of 
a,  all  combinations  of  SNR  and  D  above  the  curve  (higher  SNR 
or  lower  D)  result  in  acceptable  phase  lock  while  those  below 
the  curve  do  not. 

It  is  interesting  that  for  small  frequency  deviations  lock¬ 
on  occurs  more  readily  for  large  values  of  a  (rapid  frequency 
fluctuations)  than  for  small  a,  while  for  large  frequency  deviation 
signals  with  small  values  of  a  (long  slow  fluctuations)  can  be 
tracked  at  a  lower  SNR  than  can  those  with  large  a.  There  is 
good  physical  Justification  for  this.  If  the  modulation  index 


*10 


+15 
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is  small  compared  to  1,  the  system  does  not  have  to  track 
the  short  term  phase  fluctuations  and  performs  best  when  these 
fluctuations  are  rapid  and  soon  average  to  zero.  For  I 
greater  than  1,  the  system  must  follow  the  phase  changes  to 
remain  locked.  It  is  then  an  advantage  if  they  occur  slowly. 
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AN  EXAMPLE 

Since  the  natural  variables  required  to  describe  the  frequency 
estimator  performance  are  somewhat  different  from  those 
commonly  used  in  signal  processing  discussion,  an  example  of 
the  application  of  these  results  is  appropriate.  Suppose 
the  sine  wave  to  be  tracked  is  known  to  have  an  rms  frequency 
deviation  of  five  Hertz  about  its  center  frequency  with 
a  "typical”  relaxation  time  of  one  minute  for  the  fluctuations. 
Suppose  rms  amplitude  fluctuations  of  30%  are  experienced 
about  the  average  received  amplitude  and  that  these  fluctuations 
are  also  correlated  over  about  one  minute.  Then 

a  *  8  *  1/60  radian/second 

h  m  °-3 

D  ■  5  Hertz  ■  10ir  radians/second 
I  -  D/a  -  600tt  Z  1900 

to 

Notice  that  I  tends  to  be  a  rather  large  number  whenever 
frequency  fluctuations  of  this  magnitude  occur  over  periods 
of  seconds  or  minutes.  In  fact,  anytime  one  is  inclined  to 
describe  the  signal  behavior  as  a  sine  wave  of  varying  frequency 
rather  than  as  a  random  narrow  band  process,  this  pretty  much 
implies  a  value  of  I  much  larger  than  unity.  The  amplitude  modula¬ 
tion  index  IA  on  the  other  hand  is  essentially  constrained  to  be 
less  than  unity. 

It  is  usual  in  signal  processing  problems  of  this  class  to 
represent  the  signal  to  noise  ratio  (SNR)  as  the  total  signal  power 

43 


jjZ  -  -  -  -  iirf  } 


NSWC/WOL/TR  75-74 


2 

(A  /2)  divided  by  the  noise  power  in  a  one  Hertz  band,  considering 

only  positive  real  frequencies.  S^nce  the  noise  power  density  N 

in  this  treatment  was  defined  for  a  two-sided  spectrum,  the  noise 

in  a  one  Hertz  band  of  a  one-sided  spectrum  is  2N.  Thus  the 

SNR  is  equal  to  A  /4N  or  p/2.  The  nondimensionalized  signal 

to  noise  ratios  p.  and  p  are  thus 

KA  to 

PA  *  p/8  -  2(SNR)/B  *  2 (SNR)/(l/60)  =  120(SNR) 
p  *  P/a  *  120 (SNR) 

With  these  quantities  thus  defined  we  can  use  equations 
derived  in  the  previous  sections  to  determine  system  parameters 
and  performar.  le  as  a  function  of  the  SNR.  In  particular  equation 
(40)  gives  the  time  constant  of  the  amplitude  estimator  as 

tA  ■  1/b/i  +  2l/pA' 

-  60//  +  2(0.3)2(120)  (SNR) 

-  60//l  +  21. 6 (SNR)'  seconds  (56  ) 


equation  (4l)  gives  the  weigh u  for  the  measured  amplitude  information 
as 

aA  -  c/l  +  2IA2pA  -  1]//  +  2IA20A 

-  [/"  +  21. 6 (SNR)  -  l]//l  +  21. 6 (SNR )'  (57  ) 
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and  from  equation  (36)  we  obtain  the  rms  relative  error  in 
the  estimated  signal  amplitude  as 


°A/A 


{[/l  +  2IA2pA  -  13/Pa)1/2 


-  1]/120(SNR)}1/2 


(58) 


These  three  results  are  plotted  in  Figure  9  as  a  function 

of  the  SNR,  expressed  in  decibels.  Note  that  a  SNR  value 

of  about  -13  db  (that  is  SNR  -  1/21.6)  is  the  point 

at  which  the  amplitude  estimator  begins  to  work.  Below 

this  point  the  input  data  is  nearly  ignored  and  the  estimation  error 

is  nearly  equal  to  the  input  fluctuation  amplitude.  For  higher 

values  of  SNR  the  input  measurements  are  more  heavily  weighted  and 

the  time  constant  becomes  shorter  to  track  the  fluctuations 

more  precisely. 

Similarly  for  the  frequency  and  phase  estimation  part  of 
the  system,  equation  06)  gives  the  bandwidth  of  the  phase  lock 
loop  as 


the  rms  phase  error  can  be  determined  from  equation  (37)  as 


RELATIVE  AMPLITUDE  ERROR  (  cr^A A)  INPUT  DATA  WEIGHT  (GA>  TIME  CONSTANT  (r)  IN  SECONDS 


-10  0  10  20  30 

SIGNAL  TO  NOISE  RATIO  (SNR)  IN  DE< 
FIG.  9A  VARIATION  OF  TIME  CONSTANT  WITH  SI 


-10  0  10  20  30 

SIGNAL  TO  NOISE  RATIO  (SNR)  IN  DE< 

FIG.  9B  VARIATION  OF  INPUT  DATA  WEIGHT  WITH 


-10  0  10  20  30 

SIGNAL  TO  NOISE  RATIO  (SNR)  IN  DE« 

FIG.  9C  VARIATION  OF  RMS  AMPLITUDE  ERROR  V 


FIG.  9  PARAMETERS  AND  PERFORMANCE  OF  EXAMPLE  AMPLITUDE 
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«e ’ -  i]/pu)1/2 

-  {[/l  +  2(1900)/2(120)(SNR)V  -  1]/(120) (SNR) }1/2 

-  {[/l  +  (5.9  10^ )  (SNR)1^2  '  -  1]/(120)  (SNR) } 1/2  (60) 


and  the  rms  frequency  error  is  obtained  from  equation  (38)  as 

%  m  /PUU  •  +  2IU/2^'-  «-(lt  2Iu^)1/V(2pM)1/2 

-  (l/60)[/  +  (5.9  104)(SNR)1/2  -  !]• 


(1  +  5.9  101<)(SNR)1/2)1/4/(240(SNR))1/2 


«  +  (5.9  10i*)(SNR)1/2'  -  1] » (1  +  (5.9  104)  ( SNR ) 1/2 ) 1/11 

930(SNR)1/2 


(61) 


These  three  functions  are  plotted  in  Figure  10  as  a  function 
of  the  SNR,  expressed  in  decibels.  Note  that  the  requirement 


*  120(SNR) >  1/8I2  «  1/8 (1900)2 

is  satisfied  at  an  SNR  of  2.88  10”10  or  about  -96  db,  so 
for  all  reasonable  values  of  SNR  the  asymptotic  forms  in 
equations  (50)  and  (51)  hold  to  give 


°6  *  (Iw)1/2*(2/pw)3/8  “  (1900)1/2,(60(SNR))"3/8 

-  1.J|2/(SNR)3/8  radians  <62> 


J,.y,  2. ... „  .,>4 


±-3uiL 


RMS  FREQUENCY  RROR  (<ru)  IN  RAD/SEC  RMS  PHASE  ERROR  (aQ  )  IN  RADIANS  LOOP  BANDWIDTH  IN  RADIANS/SECOND 
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FIG.  10A  VARIATION  OF  LOOP  BANDWIDTH  WITH  SNR 
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FIG.  10B  VARIATION  OF  RMS  PHASE  ERROR  WITH  SNR 


FIG.  10C  VARIATION  OF  RMS  FREQUENCY  ERROR  WITH  SNR 


FIG.  10  PARAMETERS  AND  PERFORMANCE  OF  EXAMPLE  FREQUENCY  ESTIMATOR  VS.  SNR 
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ou  -  a(2)1/2(I1,)3/*(2/pu)1/8 

-  (l/60)(l.lllll)(1900)3/1,(60(SNR))“1/8 

■  4.06(SNR)1/^8  radians/second  (63) 

The  curve  for  Cg  shows  the  rms  phase  error  becoming  less  than 
one  radian  for  SNR  values  above  2.52  or  about  +4  db.  This 
can  be  considered  the  threshold  at  which  phase  lock  occurs. 
Below  this  SNR,  the  expressions  for  frequency  and  phase 
errors  are  probably  meaningless  and  are  shown  dashed  in 
Figure  10.  This  minimum  SNR  point  could  have  been  determined 
directly  from  Figure  8  using  D  *  5  Hertz  and  a  *  ,016.  Note 
that  because  of  the  very  weak  dependence  of  on  the  SNR, 
the  rms  frequency  estimation  error  is  between  about  1  and 
4  radians  per  second  from  this  lock  on  SNR  to  about  +50  db. 
Again  it  should  be  remembered  that  this  weak  dependence  on 
SNR  is  due  in  part  to  the  model  assumed  for  the  signal 
fluctuations.  It  should  also  be  remembered  that  the  curves 
in  Figures  9  and  10  are  not  for  a  fixed  parameter  system 
operating  over  a  range  of  SNR  values.  They  represent  a  system 
that  for  each  SNR  has  the  optimum  parameters  for  that  SNR 
value,  due  to  the  amplitude  dependence  of  the  gain  coefficients 
shown  in  Figure  2. 
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SUMMARY  AND  CONCLUSIONS 

This  paper  has  addressed  the  problem  of  finding  the 
optimum  frequency  estimator  for  a  sinusoid  whose  frequency 
is  fluctuating  randomly  in  time  and  whose  amplitude  is 
fluctuating;  in  the  presence  of  additive  random  noise.  This 
is  essentially  the  problem  of  demodulation  of  a  frequency 
modulated  signal  through  a  noisy,  fading  transmission  channel. 
The  desired  signal  properties  were  modeled  using  continuous 
state  variable  methods  and  the  optimum  estimator  derived  using 
a  linearized  form  of  continuous  Kalman  filter  theory.  The 
problem  is  sufficiently  small  to  permit  expansion  of  the 
matrix  equations  for  the  estimator,  to  allow  determination 
of  the  form  of  the  optimum  estimator.  Once  the  unimportant 
high-frequency  terms  are  removed,  the  resultant  system  is  shown 
to  resemble  a  phase  lock  loop  for  frequency  tracking  with  a 
coherent  detector  and  low-pass  filter  for  amplitude  estimation. 
A  mora  complex  model,  using  Rayleigh  fading  in  the  transmission 
channel,  is  covered  in  the  Appendix.  It  is  shown  that  the 
optimum  estimator  is  again  a  phase  lock  loop,  but  both  phases 
of  the  signal  amplitude  must  be  detected  and  used  to  phase- 
correct  the  phase  lock  loop  reference  signal.  In  both  cases 
it  is  satisfying  to  find  that  the  mathematically  optimum 
processor  is  essentially  the  same  as  one  would  have  been 
tempted  to  design  on  an  ad  hoc  basis. 
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It  is  also  possible  to  solve  the  differential  equations 
for  the  uncertainty  matrix  for  the  steady  state  case,  in  terms 
of  the  input  signal  and  noise  parameters.  This  permits  one  to 
determine  the  optimum  parameters,  as  well  as  the  form,  of  the 
processor  and  to  determine  the  steady  state  estimation  errors. 
Thus  for  the  steady  state  case  we  have  fully  determined  the 
design  of  the  optimum  processor  and  its  performance,  as  a 
function  of  the  characteristics  of  its  inputs.  These  can  either 
be  used  directly  in  design  of  an  estimator  or  can  serve  as  a 
reference  against  which  an  ad  hoc  processor  can  be  compared. 

The  optimum  parameters  of  the  processor  depend  upon  the 
modulation  characteristics  of  the  signal  (deviation  and  bandwidth 
of  the  frequency  modulation,  percentage  modulation  and  bandwidth 
of  the  amplitude  modulation)  and  on  the  signal  to  noise  ratio. 
Thus  any  particular  set  of  parameters  chosen  represents  a  point 
design  for  a  given  set  of  signal  characteristics.  However,  since 
the  signal  characteristics  can  all  be  measured  by  the  processor, 
the  information  derived  here  can  also  be  used  to  design  an 
adaptive  estimator  which  measures  the  characteristics  of  the 
received  signal  and  adjusts  its  parameters  accordingly. 
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APPENDIX 

EFFECT  OF  RAYLEIGH  FADING  SIGNAL  MODEL 

The  signal  model  used  in  the  main  body  of  this  report 
assumed  a  phase-continuous  signal  of  varying  frequency,  with 
random  modulation  in  amplitude.  It  was  shown  that  the  optimum 
frequency  estimator  for  this  problem  consisted  of  a  phase  lock 
loop  system.  It  can  be  argued  that  this  result  is  sensitive 
to  the  signal  model,  and  that  the  phase  lock  loop  is  tracking 
only  the  "carrier"  portion  of  the  amplitude  modulated  signal. 
Thus  the  amplitude  fluctuations  have  no  effect  on  the  resultant 
processor.  The  sorts  of  physical  phenomena  for  which  this 
signal  model  might  be  appropriate  are  those  in  which  the  signal 
source  itself  fluctuates  in  amplitude,  its  coupling  to  the 
transmission  medium  is  highly  time  or  aspect  dependent  as  with 
complex  radiation  patterns,  or  the  transmission  path  itself 
has  a  time  dependent  gain  such  as  with  moving  absorptive  cloud 
masses . 

In  another  broad  class  of  problems,  the  amplitude 
fluctuation  is  due  to  variations  in  the  transmission  path 
due  to  multipath  propagation.  In  both  sky-wave  high  frequency 
radio  propagation  and  in  underwater  acoustic  propagation,  it 
is  common  for  the  signal  to  arrive  by  several  paths  of 
approximately  equal  strength  but  whose  path  lengths  differ 
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by  some  number  of  wavelengths.  As  these  signals  Interfere  with 
each  other  the  resultant  amplitude  varies,  but  the  phase  also 
varies  due  to  the  vector  addition  of  signal  arrivals  of  randomly 
varying  phase  relationship.  The  most  common  analytical  form 
used  for  such  signals  assumes  that  each  arrival  path  of  the 
signal  consists  of  one  random  component  of  sine  phase  and  one 
of  cosine  phase.  Then  if  many  such  signals  are  added,  the 
Central  Limit  Theorem  indicates  that  the  amplitudes  of  the 
resultant  sine  and  cosine  components  are  independent  zero-mean 
Gaussian  variables.  Prom  this  model,  the  total  amplitude  may  be 
shown  to  follow  a  so-called  Rayleigh  distribution  while  the 
resultant  phase  angle  is  a  uniform  random  variable.  This 
Rayleigh  model  agrees  well  with  experimental  evidence  and 
includes  such  complications  as  rapid  reversal  of  signal  phase 
at  signal  dropouts.  Since  the  phase  is  not  continuous  and  there 
is  no  "carrier"  to  track,  one  would  expect  that  this  signal 
model  would  have  some  effect  on  the  form  of  the  optimum  receiver. 

Following  the  same  procedure  as  in  the  main  body  of  the 
report,  the  received  signal  state  can  be  represented  by  the  vector 


6 

X]  ■  u> 

A 
B 


(64) 
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such  that  the  measurable  signal  at  any  Instant  Is 

h(X)  -  h(0  ,u>,  A,B)  -  A  cos0  +  B  sin6  (65) 

A  and  B  are  thus  the  random  amplitudes  of  the  two  components 
of  the  received  signal  and  are  zero-mean  processes  of  equal 
variance.  Again  the  frequency  parameter  o>  does  not  directly 
appear  in  the  measurable  quantity. 

There  is  another  more  serious  problem  with  this  formulation* 
Since  h(X),  and  its  evolution  with  time,  can  be  represented  by 
many  combinations  of  A,  B,  and  0  (for  example  adding  it  to  the 
initial  value  of  0  and  reversing  the  signs  of  A  and  B),  the 
true  state  of  the  system  is  not  observable.  This  inherent 
coupling  of  the  uncertainties  of  the  amplitude  and  phase 
variables  eventually  causes  analytical  difficulties  in  the 
solution  of  this  problem.  However  it  is  still  possible  to 
determine  the  form  of  the  optimum  receiver. 

Again  we  need  to  represent  the  evolution  of  the  state 
vector  X  by  an  equation  of  the  form  X  ■  FX+Gw,  and  the  appro¬ 
priate  matrices  are 

F  - 
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and 


0 


0 


0 


G 


o  0  0 

0  8  0 


0 


0 


(67) 


The  w  vector  is  a  three  component  random  process  with  mean 
value  and  covariance  given  by 


w 


w. 


0 


(68) 


(69) 


These  definitions  of  F,  G,  and  w  again  establish  a  system  with 
a  source  center  frequency  of  with  random  low-pass  frequency 
modulation  of  "bandwidth"  a.  The  rms  deviation  from  the  center 
frequency  is  again  given  by  ^aQ^/2 .  The  complex  amplitude 
modulation  is  a  low-pass  Rayleigh  process  of  bandwidth  8, 
and  the  mean  power  in  the  signal  is  determined  by  the  value  of 
qa«  Note  the  mean  amplitude  and  the  rms  fluctuations  in  amplitude 
are  not  separately  specified  as  before,  since  they  are  inherently 
related  by  the  Rayleigh  distribution. 
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Additional  matrices  to  be  defined  as  Inputs  to  the  Kalman 
filter  solution  are  the  measurement  noise  matrix  R,  which  is  a 
lxl  matrix  equal  to  the  additive  noise  power  density  N,  the 
linearized  measurement  matrix  H,  which  is 


[h]«  [sh/30,  3h/3w,  3h/3A,  3h/3B  ] 

tA  A  /\  A  A  A 1 

-Asin0+Bcos0 ,  0,  cos0,  sin0 

and  the  estimation  uncertainty  matrix  P,  which  has  the  form 


(70) 


00 

P0O) 

P0A 

P0B 

» 

U)0 

P0)0> 

P0)A 

P0)B 

A0 

PAu) 

PAA 

PAB 

B0 

PBo) 

PBA 

PBB 

(71) 


From  these  matrix  definitions  it  is  now  possible  in  principle 
to  derive  the  optimum  Kalman  filter  estimator  and  to  determine 
the  performance  from  the  steady  state  values  of  the  P  matrix. 
Again  this  is  considerably  harder  than  it  sounds. 

Following  the  same  procedure  as  in  the  main  body  of  this 
report,  the  first  step  is  to  expand  the  differential  equation 
for  P,  so  that 
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P  -  P  P  +  P  FT  +  G  Q  GT  -  P  HT  R"1  H  P 


?P 

*00) 

Pa)U)”oP0w 

PuA"8P0A 

PuB’6P6B 

Pw<i>"*oP0u) 

a2Q  -2aP 

ww 

"aPwA~8Pu)A 

-«PuB-BP( 

PojA“BP0A 

"aPa)A"BPwA 

•V2®** 

-2bpab 

PwB“8P0B 

"aPa>B”8Pu)B 

-2«pab 

62Qa-2BP 

— ( 1/N ) 


P00(Bcos0-Asin0 )  +  PQAcos0  +  P0gSin0 
P0w(Bcos0-Asin0 )  +  Puccis©  +  P^gSin© 
P0A(Bcos0-Asin0)  +  PAAcos0  +  PAgSin0 
Pfln(Bcos0-Asin0)  +  PAtJcos0  v  PRnsin0 


P00(Bcos0-Asin0)  +  P0Acos0  +  P0gSin0 
P0w(Bcos0-Asin0 )  +  PuAcos0  +  P^gSinO 
P0A(Bcos0-Asin0)  +  PAAcos0  +  PAgSin0 
P0g(Bcos0-AsS.n0 )  +  PAgCOS0  *  PggSin0 


(72) 


In  the  above,  advantage  has  been  taken  of  the  fact  that  PAj  “ 
to  simplify  the  equations  slightly. 

We  next  wish  to  take  advantage  of  the  fact  that,  once  the 
indicated  multiplication  is  performed  in  equation  (72)  above, 
many  of  the  expressions  will  contain  terms  oscillating  at  2ui  . 
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Since  these  do  not  influence  the  long-term  evolution  of  the  P 
matrix,  they  will  be  dropped.  The  resultant  ten  scalar  equations 
for  the  elements  of  P  are: 
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-«PUB-6,,UB-(1/2',>  |(A  +B  )Peo,PeB+B'Pea1[,AB+PuAt>9B)-A<PUBP9B+P0UPBB)  +PUAPAB  +PUBPBB  (731) 
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Recall  that  in  the  system  studied  in  the  main  body  of 
this  report,  the  differential  equations  for  the  covariance 
terms  between  the  amplitude  and  the  frequency  and  phase 
functions  were  linear  first  order  equations  whose  solutions 
would  tend  toward  zero.  Thus  in  the  steady  state,  those  two 
covariance  terms  could  be  assumed  to  be  zero,  considerably 
simplifying  the  solution.  We  shall  investigate  whether  this 
simplification  applies  to  the  present  problem. 

First  consider  equation  ( 7 3 J )  for  PAB<  If  we  collect  all 
the  terms  containing  PAB  on  the  right  hand  side,  we  obtain 

<^AB>  '  >PAB-(1/2N>  [(A2+B2)P9AP9B+BPAftP6B-APeAPBB]  (74) 

This  is  again  a  linear  first  order  differential  equation  in  PAB, 
with  an  input  driving  term  containing  A,  B,  PQA,  P0B,  PAA  and  PBB. 
While  its  homogeneous  solution  decays  to  zero,  its  complementary 
solution  depends  on  the  driving  inputs. 

Taking  a  similar  look  at  the  pair  of  equations  (73f)  and 
(73h)  for  PQA  and  PwA,  they  can  be  written 

<P6A>  -  '  )P8A+PuA-<1/2N>[B(P9ePAA+PeA2>-APe9PAB+PeBPAB] 

and 

<PwA>  -  '< '  '  '  >PuA-( '  ' '  >P9A-(1/2N>[BPe/AA-APe.PAB+PUBPAB] 

Again  this  pair  of  coupled  equations  has  a  homogeneous  solution 
which  decays  to  zero,  but  also  may  have  a  complementary  solution 


(75a) 

(75b) 
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driven  by  some  other  terms.  A  similar  situation  occurs  for  the 
two  equations  for  P0B  and  Pw0. 

On  inspection  of  this  set  of  five  equations,  it  is  seen  that 
even  if  all  five  covariances  PAB,  P0A,  P0B,  PwA,  and  PwB  were 
initially  zero,  there  are  driving  terms  in  their  differential  equations 
which  will  excite  some  non-zero  response.  An  example  of  such  a 
term  is  the  (B/2N) (P00PAA+PQA2)  term  in  equation  (75a).  Since 
the  second  factor  is  clearly  positive  and  non-zero,  the  covariance 
element  PQA  obviously  will  become  non-zero  whenever  B  (the  sine 
component  of  the  input  signal)  is  non-zero.  This  is  one  unfortunate 
result  of  the  non-observability  of  the  state  of  the  system. 

Physically,  whenever  B  is  non-zero,  there  is  no  way  to  distinguish 
between  a  change  in  the  other  signal  component  amplitude  A  or  a 
true  change  in  the  phase  0.  Thus  the  covariance  between  A  and  6 
will  be  non-zero. 

These  five  covariance  terms  are  different  from  the  remainder 
of  the  covariance  matrix  though,  in  that  they  all  are  driven  by 
A  and  B  which  are  zero-mean  random  variables.  Thus  while  the 
covariance  values  are  not  actually  zero,  they  only  fluctuate 
about  zero  in  response  to  the  A  and  B  state  inputs.  Thus  over 
the  long  term,  their  steady  state  value  averages  to  zero.  Inclusion 
of  these  terms  unquestionably  adds  to  the  complexity  of  the  state 
estimator  and,  since  they  are  dynamically  changing  parameters,  forces 
the  implementation  of  the  differential  equations  for  at  least  this 
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part  of  the  P  matrix.  However  since  they  are  due  largely  to 
the  over-specification  of  the  state  of  the  system,  it  is  felt 
that  they  can  be  dropped  from  at  least  a  first  order  approximation 
to  the  optimum  estimator.  The  estimator  can  then  be  designed 
using  only  those  P  matrix  coefficients  which  have  a  non-zero  steady 
state  value  and  do  not  require  continuous  computation. 

On  the  basis  of  this  assumption,  the  steady  state  values 
for  the  remaining  terms  of  the  P  matrix  may  be  found  by  dropping 
the  five  P  matrix  terms  discussed  above  and  setting  the  remainder 
of  the  derivitives  to  zero.  This  leads  to  the  set  of  four 
equations 


<*69* 

*  2Peu,-(A2+B2)pee2/2N  *  0 

(76a) 

<P  > 
0)0) 

“  “V2aP*-(A2+B2)peu2/2N  ■  0 

(76b) 

<pfl  > 
00) 

■  p«-oP9U-(A2+B2)peepeu/2N  "  0 

(76c) 

A 

T)  • 
> 

> 

V 

,  <PBB>  -  B2Qa-26Paa-Paa2/2N  -  0 

( 7  6d ) 

where  it  is  recognized  that  the  equations  for  PAA  and  P0B  are 
identical.  Comparing  these  equations  with  equations  (19)  in  the 
main  text  shows  that  they  are  identical  in  form,  except  for  the 
replacement  of  A  /2N  with  (A  +B  )/2N  for  the  signal  to  noise  ratio. 
Thus  the  solutions  for  the  steady  state  frequency  and  phase 
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estimation  errors  obtained  In  the  main  text,  and  the  conclusions 
reached  regarding  performance  as  a  function  of  signal  to  noise 
ratio  and  other  parameters,  remain  unchanged.  The  Interpretation 
of  the  solution  for  the  amplitude  estimation  error  must  be 
changed  slightly  because  QA  Is  now  related  to  the  total  signal 
power  rather  than  to  an  amplitude  modulation. 

The  next  step  in  finding  the  form  of  the  optimum  estimator 
is  to  substitute  the  matrix  definitions  into  the  estimation 
equation 


X  -  P  X  +  Q  wQ  +  P  HT  R"1  (z  -  h(X) ) 


0) 

A 

-au) 

A 

-BA 

A 

-BB 


0 

a«o{ 

0 

0 


+  (1/N)(z-Acos0-Bsin0)* 


P00(Bcos0-Asin0)  +  P0Acos0  +  P0Bsin0 
P0w(Bcos0-Asin0)  +  PuAcos0  +  PuBsin0 


(77) 


PQA(Bcos0-Asin0)  +  PAAcos0  +  PA0sin0 

P0B(Bcos0-Asin0)  +  PABcos0  +  PBBsin0 
■  m 

This  is  the  set  of  differential  equations  which  must  be  solved 
if  all  P  matrix  values  are  assumed  to  be  non-zero.  However,  if 
we  retain  only  the  five  terms  P00,  P uhJ,  P0(jJ,  PAA,  and  PBB  which 
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have  non-zero  steady  state  values,  and  ignore  the  remaining 
covariance  terms  which  fluctuate  around  zero,  these  equations 


simplify  to 


+ (1/N) (z-Acos0-Bsin6 ) • 


P00 (Bcos0-Asin0 ) 
PQ  (Bcos0-Asin0 ) 

DO) 

PAAC0S6 

PBBSin8 


(78) 


Finally,  if  we  again  take  short  time  averages  to  eliminate  the 
effects  of  terms  at  the  carrier  frequency  or  higher,  we  obtain 

A 

•  A  A  A  A  A 

<0>  *  a)  +  (PQ0/N)<(Bcos0  -  Asin0)*z> 


(79a) 


<d)>  *  a(w  -w)  +  (P.  /N)<(Bcos0  -  Asin0)*z> 
o  yoj 


<A> 


-6A  +  (pAA/N)  (<z-cos0>-A/2) 


(79b) 

(790 


<B>  « -BB  +  (Pgg/N)  (<z  •  sin0>-B/2 ) 


(79d) 


Figure  11  shows  this  set  of  equations  expressed  in  a  block 
diagram  form  similar  to  that  of  Figure  1  in  the  main  text .  The 
upper  part  of  the  figure  shows  the  estimators  for  the  amplitude 
components  A  and  B.  Each  is  a  simple  low-pass  filter  averaging 
the  mean  product  of  the  input  signal  with  either  the  sine  or  the 
cosine  phase  output  of  the  reference  function  generator.  The  time 
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FIG.  1 1  FREQUENCY  TRACKER  STRUCTURE  FOR  RAYLEIGH  FADING  CASE 
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constant  of  this  filter  is  determined  by  the  value  of  the 
uncertainty  matrix  element  PAA  or  PgB,  just  as  in  the  system 
studied  in  the  main  text.  In  fact  the  only  change  due  to  the 
Rayleigh  fading  model  is  that  two  channels  of  amplitude  estimator 
are  required  rather  than  just  one. 

The  frequency  and  phase  estimation  portion  of  the  circuit 
is  also  quite  like  that  of  Figure  1  in  the  main  text,  in  that 
it  is  a  second-order  phase  lock  loop  whose  dynamics  are  determined 
by  the  coefficients  P0Q  and  PQw.  Since  these  are  determined  by  the 
same  equations  in  both  of  the  models,  the  performance  of  the  phase 
lock  loop  should  be  identical. 

The  only  real  difference  caused  by  the  Rayleigh  fading  model 
is  in  the  reference  function  used  by  the  phase  lock  loop.  In  the 
simpler  case,  the  PLL  reference  came  directly  from  one  phase 
of  the  function  generator  output.  In  the  Rayleigh  case,  the 
function  Bcos0  -  Asin0  is  used  to  form  the  reference  signal  for 
the  phase  lock  loop.  Whenever  A  and  B  (or  their  estimates)  are 
constant,  this  difference  is  immaterial  and  the  Rayleigh  solution 
degenerates  to  that  of  the  previous  case.  However,  when  A  and 
B  fluctuate  as  assumed  in  the  Rayleigh  model  and  these  fluctuations 
are  tracked  by  the  A  and  B  estimators,  the  complex  reference 
function  automatically  changes  its  phase  to  compensate  for  the 
new  ratio  of  A  and  B.  Thus  phase  changes  caused  by  fading  of 
the  signal  phase  components  are  handled  by  the  amplitude  estimation 
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network  and  do  not  require  a  response  from  the  phase  and 
frequency  estimation  circuit.  This  is  particularly  important 
in  the  case  of  deep  signal  fades  which  cause  a  rapid  reversal  in 
the  received  signal  phase,  since  it  is  these  rapid  phase  changes 
which  are  the  most  difficult  for  the  phase  lock  loop  to  handle. 

Obviously  things  are  not  quite  this  simple,  again  because 
of  the  non-observability  of  the  state.  Whenever  a  change  is 
observed  in  the  received  signal  phase,  it  can  be  interpreted 
as  either  a  frequency  change  or  as  a  change  in  the  signal 
vector  components.  The  degree  to  which  this  is  handled  by 
the  amplitude  or  the  phase  portion  of  the  estimator  system  is 
determined  by  the  values  a,  6,  and  used  in  modeling  the  signal 
behavior.  The  fact  that  both  portions  of  the  system  attempt  to 
respond  to  each  change  is  the  source  of  the  fluctuating  covariance 
terms  such  as  PQA  which  we  chose  to  ignore  in  developing  this 
estimator. 

The  conclusion  to  be  reached  from  this  Appendix  is  that, 
at  least  to  a  first  approximation,  the  optimum  frequency  estimator 
in  the  presence  of  Rayleigh  fading  is  still  a  phase  lock  loop 
system.  It  is  somewhat  more  complex  than  in  the  pure  amplitude 
fading  model  in  that  two  amplitude  estimators  are  formed  from 
quadrature  phases  of  the  reference  signal,  and  these  estimators 
are  used  to  form  the  reference  signal  for  the  phase  lock  loop. 

The  primary  effect  of  this  change  is  that  it  assists  the  phase 
lock  loop  in  tracking  the  signal  phase  through  deep  fading  periods. 
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GLOSSARY 


A  signal  amplitude 

Aq  mean  amplitude  of  signal 

a  bandwidth  of  frequency  fluctuations 

B  bandwidth  of  amplitude  fluctuations 

B  amplitude  of  quadrature  signal  component 

BW  bandwidth  of  the  phase  lock  loop 

D  rms  frequency  deviation  in  radians/sec 

6(t)  Dirac  delta  function 

F  matrix  describing  evolution  of  signal  state 

G  matrix  coupling  random  process  to  signal  state 

Ga  input  signal  weight  in  amplitude  estimator 

H  linearized  measurement  matrix 

h(X)  scalar  measure  of  signal  state 

IA  amplitude  modulation  index 

I  frequency  modulation  index  «  D/o 

K  a  parameter  of  the  phase  lock  loop 

N  noise  power  per  Hertz  (double  sided) 

P  covariance  matrix  of  estimation  error 

Q  covariance  matrix  of  W 

R  noise  autocorrelation  matrix 

p  a  signal  to  noise  density  ratio  parameter 

pA  normalized  signal  to  noise  ratio  ■  p/B 

normalized  signal  to  noise  ratio  ■  p/o 
SNR  signal  to  noise  ratio  in  1  Hertz  band  (single  sided) 

T(s)  transfer  function  of  phase  lock  loop 
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X 
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GLOSSARY  (Cont.) 

time  constant  of  amplitude  estimator 
signal  phase  angle  in  radians 
additive  noise  process 

random  process  producing  fluctuations  in  state 
mean  value  of  W 

signal  frequency  in  radians/sec 
mean  frequency  of  signal 
signal  state  vector 
derivative  of  signal  stats  vector 
estimate  of  signal  state  vector 
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